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ABSTRACT 


A  series  of  spherical  one-dimensional  computer  code 
calculations  was  performed  to  assess  the  effects  of  mechani¬ 
cal  material  property  parameters  on  wave  propagation  in  dry 

/ 

porous  and  compacted  tuff.  Peak  stress  attenuation  was  found 
to  be  most  sensitive  to  initial  porosity.  In  compacted  material 
peak  stress  attenuation  was  found  to  be  more  sensitive  to  the 
shear  modulus  than  to  the  yield  strength.  Increasing  hoop 
tensile  stresses  occurred  with  increasing  yield  strength. 
Introduction  of  a  hoop  tensile  fracture  criterion  did  not 
affect  peak  stress  attenuation  but  did  affect  the  radial  dis¬ 
placement.  Generally  radial  displacements  of  particles  at 
the  same  initial  radius  from  the  source  increased  with  in¬ 
creasing  porosity.  Shock  wave  arrival  times  at  various 
distances  from  the  source  increased  with  increasing  porosity. 


I.  INTRODUCTION 


In  order  to  satisfactorily  predict  stress  wave  propa¬ 
gation  in  the  earth  resulting  from  an  underground  explosion, 
it  is  necessary  to  know  certain  material  properties  of  the 
earth  which  govern  its  response  to  dynamic  loads.  The  present 
study  was  undertaken  as  a  preliminary  attempt  to  determine 
the  relative  importance  of  various  mechanical  material  para¬ 
meters  on  spherically  symmetric  wave  propagation  in  a  rock 
which  may  be  porous.  The  problem  studied  consists  of  a 
spherical  cavity  in  the  rock  with  an  initial  radius  of  5  cm. 

At  time  t  =  0 ,  a  normal  stress  a(rc,t)  is  applied  to  the 
interior  surface  of  the  cavity.  A  one-dimensional  Lagrangian 
finite-difference  computer  code  is  used  to  compute  the 
subsequent  stress  field  in  the  rock.  The  radial  stress, 

°r(r»  *0 ,  the  tangential  stress,  a&(r,  t) ,  and  the  hydrostatic 
pressure,  P(r,  t) ,  are  calculated  as  functions  of  time  along 
the  paths  of  selected  particles  which  are  initially  at  7.5,  10, 
12,5,  25  and  30  cm  from  the  center  of  the  cavity.  Although 
all  problems  were  run  using  the  above  dimensions,  scaling 
of  the  problems  to  different  dimensions  would  not  alter  the 
relative  importance  of  wave  propagation  on  the  various 
material  property  parameters  considered.  In  addition  to  the 
stress  field,  particle  displacements  are  also  computed. 

The  mechanical  material  property  parameters  which  are 
varied  in  this  parameter  study  were  chosen  to  be  representa¬ 
tive  of  dry  Nevada  Test  Site  (NTS)  tuff.  Inasmuch  as  this 
was  a  parameter  study,  no  attempt  was  made  to  simulate  speci- 
ficially  in  detail  any  of  the  many  NTS  tuffs  which  exist  in 
the  field.  The  effect  of  moisture  was  not  considered  in 
this  study. 

The  following  sections  describe  the  material  property 
model,  the  computation  technique  and  the  results  obtained  from 
several  computations  in  which  the  values  of  the  mechanical 
parameters  were  varied. 
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2.  EQUATION  OF  STATE  MODEL  AND  CALCULATION  METHOD 


The  equation  of  state  model  used  to  represent  dry 

porous  tuff  in  the  present  parameter  study  is  the  S3  Porous 
f  1  21 

EOS  model.1  *  1  This  model  includes  a  hydrostatic  pressure- 
volume  -energy  (P-V-E)  relationship  which  provides  for  irrever¬ 
sible  void  collapse  plus  deviatoric  stresses  to  provide  for 
shear  strength.  The  P-V-E  model  is  quite  similar  to  that 
developed  by  Herrmann. ^ 

2.1  HYDRODYNAMIC  EQUATION  OF  STATE 

The  solid  (void  free)  material  is  assumed  to  obey  a 
Mie-Gruneisen  equation  of  state 

P  -  Pr  +  Gp(E  -  Er)  (1) 

where  Pr  is  some  reference  curve  in  the  P-V  plane,  Er  is 
the  specific  internal  energy  along  that  curve,  G  is  Gruneisen's 
ratio  which  is  assumed  constant,  and  p  =  1/V  is  the  density 
of  the  solid  material.  The  reference  curve  is  frequently  taken 
to  be  the  principal  Hugoniot  curve  of  the  solid  material  in 
the  form 


PH(p)  =  An  +  Bn2  +  Cn3 


(2) 


where  n  =  p/p  -  1  and  A,  B,  and  C  are  constants.  The 
o 

specific  internal  energy  along  the  Hugoniot  is  then  given  by 


EH  *  E,  +  7  PH  (V„  -  V) 


(3) 


Preceding  page  blank 
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When  dealing  with  a  porous  material  a  distension  ratio, 
a,  defined  as 


PS(P,  E) 

01  5  yPTE) 

is  introduced.  Here  pg  is  the  density  of  the  void  free 
material  at  a  given  value  of  pressure  and  internal  energy 
and  Pp  is  the  density  of  the  porous  material  at  the  same 
value  of  the  pressure  and  internal  energy.  The  compaction 
or  crushup  process  of  the  porous  material  is  then  modeled 
through  the  parameter  a.  In  general  a  could  depend  upon 
the  deviatoric  stresses  in  addition  to  the  hydrostatic  pressure 
and  internal  energy.  The  Mie-Gruneisen  equation  of  state  for  the 
porous  material  is  obtained  by  combining  Eqs.  (1),  (2)  and  (4) » 

P  =  ph  (otpp)  +  G  a  Pp(E  -  Eh)  (5) 

In  the  present  work  a  is  taken  to  depend  linearly 
on  the  hydrostatic  pressure  as  shown  in  Fig.  1.  The  initial 
distension  ratio  for  the  virgin  material  is  a  .  For 
moderate  compressions  in  which  the  pressure  remains  less  than 
a  critical  pressure  P0e>  the  material  behaves  reversibly, 
i.e.,  upon  release  it  will  recover  its  initial  distension. 

In  this  region 

a  =  a  +  a  P  P<  PB_  (6) 

oooi  -  oe 

If  the  pressure  PQe  is  exceeded,  irreversible  compaction 
of  the  material  will  occur  and  so  long  as  the  material  con¬ 
tinues  to  crush 


a  =  a'  +  a'P 

o  1 


(7) 
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(d)  sjnsssjd 
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until  all  the  voids  are  collapsed  at  which  time  a  =  1.0. 

Once  a  attains  a  value  of  unity  it  remains  at  that  value 
thereafter.  If  during  the  crushup  process  the  material  is 
unloaded  from  some  partially  compacted  state  ^ag,  Pg j  the 
unloading  will  occur  along  a  curve  such  that 

a  =  a  +  a  P  (8) 

0  1 

This  a-P  curve  is  reversible  in  that  any  reloading  of  the 
material  will  occur  along  this  curve  until  the  state  (ae»  ?e) 
is  once  again  reached.  Upon  reaching  (ae>  Pg)  subsequent 
loading  occurs  along  the  irreversible  curve,  given  by  Eq.  (7). 
Thus,  for  a  material  which  has  been  partially  compacted  and 
then  unloaded,  the  pressure  Pg  becomes  the  new  reversible 
limit  just  at  P0g  was  the  reversible  limit  in  the  virgin 
material.  The  reversible  and  irreversible  branches  of  the 
a-P  curves  are  denoted  by  arrowheads  in  Fig.  1. 

The  degree  of  irreversible  compaction  at  any  point  in 
the  crushup  process  is  accounted  for  by  a  progress  parameter 
<f>.  This  parameter  varies  from  one  to  zero  as  the  material 
passes  from  the  porous  state  a  =  a  e,  where  irreversible 
crushup  begins,  to  the  state  a  =  1  where  compaction  is 
complete . 


For  example,  if  a  =  2.0,  and  at  some  partially  crushed 

state  a  =  1.5,  then  <J>  =  0.5  denoting  that  the  material  is 

halfway  crushed.  The  coefficients  a  and  a  in  Eq.  (8) 

0  1 

for  the  reversible  loading-unloading  path  from  a  partially 
compacted  state  are  computed  from  the  coefficients  and 

a  in  Eq.  (6)  for  the  initial  reversible  path  and  from  the 

0  l 
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progress  parameter  <f>.  These  coefficients  are  given  by 

a  =(a  -  l)  <J>  +  1 

0  '00  1 

(10) 

“i  =  \,Ke/P(l») 

Thus  if  only  a  small  amount  of  permanent  crushing  has  occured, 
i.e.,  (a  ,  P  )  is  very  near  a  ,  P  ),  the  loading-unloading 
a-P  curve  will  be  very  much  like  that  of  the  virgin  material, 
hwereas  if  compaction  is  almost  complete  it  will  be  very  nearly 
the  vertical  line  a  =  1.  The  progress  parameter,  <J> ,  varies 
along  the  irreversible  crush  curve  and  remains  constant  on  the 
reversible  curves.  The  branches  of  the  hydrostatic  curve  in 
the  plane  corresponding  to  the  branches  of  the  a-P  curve 
are  shown  in  Fig.  2.  Pcr  is  the  pressure  at  which  a  reaches 
unity. 

2.2  MATERIAL  STRENGTH 

In  addition  to  the  hydrostatic  pressure,  which  is  com¬ 
puted  from  the  equation  of  state  model  just  described,  rigidity 
and  shear  yielding  affect  wave  propagation  in  the  porous 
material.  The  principal  stresses,  o^,  principal  stress  devia¬ 
tors,  S^,  and  the  hydrostatic  pressure  are  related  by 

oi  =  -p  +  S.,  i  “  1,  2,  3  (11) 

In  the  region  where  the  material  behaves  linearly  elastic 

-  2v(h  •  y  v)  tl2) 

where  the  dots  denote  time  derivates,  p  is  the  shear  modulus, 
and  is  a  principal  strain.  In  the  porous  EOS  model  pro¬ 

vision  is  made  to  allow  the  shear  modulus  to  vary  as  the 
material  is  being  crushed.  The  effective  shear  modulus  p 
at  any  point  during  crushup  is 


ye  =  4>y p  +  (1  -  <J>) 


(13) 
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where  pp  and  pg  are  respectively  the  shear  modulii  of 
the  porous  and  solid  materials  and  <j>  is  the  crushup  progress 
parameter  discussed  earlier. 

The  stress  deviators  are  limited  by  the  von  Mises  yield 
condition 


S*  +  S2  +  S2  <  l  Y2  (14) 

1  2  3  3 

where  Y  is  a  yield  strength  which  may  depend  upon  the  hydro¬ 
static  pressure.  In  the  present  work  its  dependence  upon 
pressure  appears  only  through  the  crushup  progress  parameter. 

As  with  the  shear  modulus,  the  yield  strength  Y  in  the  par¬ 
tially  compacted  region  is  given  by 

Y  -  4>  Y#p  +  (1  -  40  Yos  (15) 

where  YQp  and  Yqs  are  respectively  constant  values  of  the 
yield  strength  of  the  virgin  porous  and  of  the  completely  compacted 
material.  No  pressure  enhancement  of  the  yield  strength  of  the 
compacted  rock  is  considered. 

2.3  METHOD  OF  COMPUTATION 

The  spherical  wave  propagation  problems  in  porous  and 
compacted  dry  tuff  were  solved  using  a  one-dimensional  Lagrangian 
elastic-plastic  computer  code  called  WORKHORSE .  ^  This  code 
solves  the  equations  of  conservation  of  mass,  momentum  and 
energy  using  the  artificial  viscosity  technique  to  smooth 
shocks.  The  problem  consists  of  a  sphere  of  tuff  which  is 
initially  35  cm  in  radius  and  has  a  5-cm  spherical  cavity  in 
the  center.  Ten  problems  were  run  in  which  an  exponentially 
decaying  normal  stress  pulse  of  the  form 

o(r  ,  tj  =  35e’t//^  kbar  (t  in  ysec)  (10) 


9 


was  applied  to  the  interior  of  the  cavity.  Two  problems  were 
run  in  which  a  square  pulse  35  kbar  for  10  Msec  was  applied 
to  the  interior  of  the  cavity.  The  response  of  the  material 
as  a  function  of  time  was  computed  at  stations  which  were  ini¬ 
tially  7.5,  10,  12.5,  25,  and  30  cm  from  the  center  of  the 
sphere.  For  computational  purposes  the  entire  35-cm  radius 
sphere  is  divided  into  200  zones  with  the  central  cavity  as 
zone  one.  The  remaining  199  zones  are  chosen  such  that  the 
mass  of  each  zone  increases  by  a  factor  of  1.02,  i.e., 

=  1.02.  Once  the  zone  size  is  established  the  time 
step  in  the  calculation  is  controlled  by  the  Courant  condition 
which  depends  upon  the  zone  size,  the  sound  speed,  and  the 
artificial  viscosity  parameters. 

At  each  cycle  in  the  calculation,  values  of  the  hydro¬ 
static  pressure,  stress  deviators,  displacement,  velocity  and 
the  time  are  printed  for  the  zones  being  monitored.  The  stress 
field  results  are  put  on  tapes  from  which  plots  are  obtained 
for  the  hydrostatic  pressure,  radial  and  tangential  stresses 
as  functions  of  time  for  the  monitored  zones. 

The  values  of  the  equation  of  state  parameters  used 
in  the  calculations  and  the  results  of  the  calculations  are 
summarized  in  the  next  section. 


3.  SPHERICAL  WAVE  PROPAGATION  IN  TUFF 


In  all  of  the  spherical  problems  run  in  this  study, 
the  Hugoniot  curve  Pjj(V)  of  the  solid  compacted  material 
was  taken  to  be 

PH(n)  =  An  +  Bn2  n  -  -  l 

M  o 

where  A  =  245.7  kbar  and  B  =  298.7  kbar.  Although  void- 
free  dry  tuff  does  not  occur  naturally,  this  estimate  for  the 
pressure-volume  relationship  was  made  by  Hastings ^  based  on 
the  constituents  of  tuff.  The  solid  density  at  zero  pressure 
is  2.4  g/cm3.  Also  throughout  all  calculations,  the  value  of 
Gruneisen's  ratio  G  was  taken  to  be  0.33. 

The  effect  on  wave  propagation  of  changes  in  the 
following  parameters  was  investigated: 

1.  pcr»  the  hydrostatic  pressure  at  which  all 
voids  are  compacted 

2.  Y,  the  yield  strength  in  the  von  Mises  yield 
condition 

3.  y,  the  shear  modulus 

4.  a  ,  the  initial  porosity 

o  o 

5.  The  effect  of  fracture  due  to  hoop  tensile  stress 

6.  Variation  of  the  elastic  ct 1  s  in  the  a-P 
mode  1 . 

The  values  of  the  parameters  used  in  each  of  the  twelve  com¬ 
puter  runs  are  presented  in  Table  I.  In  all  except  two  cases 
the  stress  applied  to  the  interior  of  the  cavity  was  the 
exponential  decay  form  given  in  Eq.  (16).  The  two  cases  in 
which  the  square  pulse  was  applied  are  indicated  in  Table  I. 
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PARAMETERS  USED  IN  SPHERICAL  CALCULATIONS 
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3.1  PEAK  STRESS 


Peak  radial  stress,  a  ,  peak  tangential  stress,  aQ, 
and  peak  hydrostatic  pressure,  P,  as  functions  of  radius,  are 
given  in  Table  II  •  Peak  tensile  stresses  and  precursor  ampli¬ 
tudes  are  also  indicated  in  some  of  the  runs  in  Table  II. 

Curves  of  the  attenuation  of  peak  radial  stress  as  a  function 
of  radius  are  plotted  for  selected  runs  in  Fig.  3.  All  runs 
were  not  plotted  as  the  results  of  several  were  quite  similar. 
However,  the  curves  shown  in  Fig.  3  do  span  all  of  the  radial 
stress  attenuation  curves  for  the  ten  runs  in  which  the  exponen 
tial  decay  stress  pulse  was  applied  to  the  cavity.  In  all 
cases  the  attenuation  was  lowest  in  the  compacted  material, 
greater  in  the  5-percent  porous  material  and  greatest  in  the 
20-percent  porous  material.  This  variation  of  the  porosity 
from  zero  to  20  percent,  leaving  the  other  parameters  constant 
(compare  runs  0D,  1C,  and  3C) ,  gave  greater  peak  stress  attenua 
tion  than  variation  of  the  other  material  parameters  over  a 
range  which  might  reasonably  be  representative  of  dry  tuff. 

Five  runs  were  made  on  compacted  material  to  assess 
the  effect  on  attenuation  of  parameters  other  than  porosity. 

The  material  in  run  0D  (y  =  223  kbar,  Y  =  1.0  kbar)  may  be 
considered  as  a  standard  and  the  radial  stress  attenuation 
curve  for  this  material  is  shown  in  Fig.  3.  The  material  in 
run  5A  was  the  same  as  0D  except  that  Y  was  increased  to 
10  kbar.  The  rate  of  attenuation  of  peak  stress,  not  shown 
in  Fig.  3,  was  slightly  greater  than  in  run  OD  (see  Table  II. 
However,  greater  hoop  tensile  stresses  and  larger  precursors 
occurred.  Run  6A  was  identical  to  5A  except  that  the  hoop 
stress  was  set  equal  to  zero  after  attaining  a  value  of  1  kbar 
in  tension.  The  simulated  fracture,  occurring  in  tension, 
produced  no  effect  on  the  attenuation  of  peak  stress,  but 
did  naturally  limit  the  tension  occurring  after  the  peak, 
and,  as  we  shall  see,  changed  the  displacement  and  velocity 


PEAK  STRESS  AS  A  FUNCTION  OF  RADIUS 
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NEGATIVE  VALUES  DENOTE  PEAK  TENSILE  STRESS 
NUMBERS  IN  PARENTHESIS  DENOTE  PRECURSOR  AMPLITUDE 
SECOND  MAXIMUM  OCCURRING  AFTER  TENSION 


field  from  that  when  no  fracture  occurred.  In  run  8B  the 
fracture  condition  was  removed,  the  yield  Y  was  kept  at 
10  kbar  and  the  shear  modulus  y  was  reduced  to  45  kbar. 

As  can  be  seen  in  Fig.  3  decreasing  the  shear  modulus  by 
a  factor  of  five  reduces  the  attenuation  enough  to  more  than 
compensate  for  the  increase  in  attenuation  produced  by  in¬ 
creasing  the  yield  strength  by  a  factor  of  ten,  e.g.,  run  0D 
to  5A.  In  problem  9A  in  which  the  square  pulse  was  applied 
to  the  cavity  interior,  the  material  parameters  were  all  the 
same  as  those  in  problem  0D.  The  attenuation  of  the  peak  stress 
in  problem  9A  was  much  slower  than  in  0D  as  can  be  seen  from 
Table  II.  Radial  stress-time  profiles  for  problem  OD  are  pre¬ 
sented  in  Fig.  4. 

Five  problems  (four  with  exponential  decay  pulse  and 

one  with  the  square  pulse  applied  to  the  cavity  interior) 

were  run  for  five  percent  (a  =  1.05)  porosity  material. 

oo 

The  material  in  the  run  1C  (y  =  45-223  kbar,  Y  =  0. 1-1.0  kbar, 
Pcr  =  20  kbar)  may  be  considered  as  a  standard.  The  attenua¬ 
tion  of  peak  radial  stress  for  this  run  is  shown  in  Fig.  3. 
Problem  2A  was  run  with  the  same  parameters  as  1C  except  that 
the  crushup  pressure,  P  ,  was  increased  by  a  factor  of  two 
to  40  kbar.  From  Table  I  and  Fig.  3  it  is  apparent  that  in¬ 
creasing  the  crushup  pressure  reduces  the  rate  of  peak  stress 
attenuation.  In  run  4A  all  was  the  same  as  in  run  1C  except 
that  the  range  of  the  yield  strength  Y  was  increased  by  an 
order  of  magnitude.  Here,  as  in  the  compacted  material,  the 
attenuation  rate  of  peak  stress  was  increased.  In  run  7A  all 
parameters  were  the  same  as  standard  run  1C  except  that  the 
a  coefficients  for  the  elastic  portion  of  the  hydrostat  were 
set  equal  to  those  of  the  irreversible  portion.  No  significant 
change  in  peak  stress  attenuation  was  observed  from  that  of 
run  1C.  This  is  apparently  due  to  the  fact  that  at  five 
percent  porosity,  the  elastic  portion  of  the  hydrostat  and 
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the  crushup  portion  in  1C  had  very  nearly  the  same  slope  so 
that  varying  the  a  coefficients  as  was  done  in  7A  had  little 
effect  on  the  release  curves.  Radial  stress-time  profiles  at 
the  monitored  stations  for  run  1C  are  presented  in  Fig.  5. 

For  comparison  of  the  effect  of  the  shape  of  the  pulse  applied 
to  the  cavity  interior,  radial  stress-time  profiles  for  run  10A, 
which  is  identical  to  1C  except  that  a  square  pulse  is  applied, 
are  presented  in  Fig.  6. 

In  run  3C,  the  initial  porosity  is  increased  to 
ot^  =  1.25.  This  run  exhibited  the  greatest  attenuation  of 
peak  stress  with  distance  as  can  be  seen  from  the  radial  stress 
attenuation  curve  in  Fig.  3  and  the  tabulated  values  in  Table  II. 
Radial  stress-rtime  profiles  for  run  3C  are  presented  in  Fig.  7. 
Run  11A  was  the  same  as  3C  except  that  the  a  coefficients 
in  the  elastic  portion  of  the  hydrostat  were  set  equal  to 
those  of  the  crushup  portion.  In  this  case  the  attenuation 
of  peak  stress  was  slightly  less  than  that  of  3C.  This  re¬ 
sult  is  what  one  would  anticipate  since  the  release  curves  are 
more  parallel  to  the  crush  curve  in  run  11A  than  in  3C.  Having 
the  release  curves  more  parallel  to  the  crush  curves  simulates 
a  more  reversible  release,  i.e.,  more  opening  up  of  pores  on 
release,  thereby  producing  less  dissipation  of  energy  and  less 
attenuation.  In  the  five-percent  porosity  cases,  runs  1C 
and  7A,  this  effect  was  also  anticipated  but,  as  mentioned 
earlier,  the  difference  in  the  release  curves  in  the  two 
cases  apparently  was  not  significant  for  those  two  cases. 

3.2  ARRIVAL  TIMES  AND  PULSE  WIDTHS 


Arrival  times  and  pulse  widths  at  1/2  peak  and  1/4  peak 
stresses  are  presented  in  Table  III.  Some  of  the  time  of -arrival 
(TOA)  columns  contain  two,  three  or  four  entries.  The  first 
entry  always  refers  to  the  time  of  first  motion.  The  last 
entry  is  always  the  arrival  time  of  the  peak.  In  some  cases 
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only  a  ramped  precursor  was  evident  and  in  those  cases  three 
entries  appear,  the  middle  one  being  the  arrival  time  of  the 
bottom  of  the  plastic  wave.,  In  other  cases  a  more  well-defined 
precursor  was  formed  and  in  those  cases,  denoted  by  four 
entries,  the  second  entry  is  the  arrival  time  of  the  precursor 
peak  and  the  arrival  time  of  the  bottom  of  the  plastic  wave. 

As  can  be  seen  from  the  entries  in  Table  III,  the  arrival  times 
for  the  waves  at  the  different  stations  are  generally  lowest 
in  the  compacted  material,  greater  in  the  five-percent  and 
greatest  in  the  20-percent  porous  material.  The  pulse  widths 
at  1/2  peak  and  1/4  peak  show  remarkably  little  variation  from 
station  to  station  in  a  given  run  in  many  of  the  runs  although 
the  peak  stresses,  as  we  have  seen,  have  greatly  attenuated. 

3.3  DISPLACEMENT  AND  VELOCITY 


The  position  and  velocity  of  the  monitored  grid  cells 
at  about  75  ysec  after  application  of  the  pulse  to  the  cavity 
interior  are  presented  in  Table  IV.  The  problems  on  the  compacted 

material  generally  ran  in  fewer  cycles  corresponding  to  a 

shorter  total  wave  propagation  time  than  the  problems  on  the 
porous  material.  Total  times  ran  from  about  75  ysec  to  over 
100  ysec.  Therefore,  in  order  to  compare  positions  and  velo¬ 
cities  at  the  same  time,  75  ysec  was  chosen.  The  particles 
initially  located  at  7.5,  10,  and  12.5  cm  from  the  cavity 
center  generally  show  increasing  displacement  and  velocity 
with  increasing  porosity.  In  the  compacted  material  runs  5A 
and  6A, where  the  yield  strength  was  increased  to  10  kbar, 
show  an  effect  of  pulling  back  toward  the  center  of  the 

problem.  In  fact,  in  run  6A,  which  differs  from  5A  only  in 

that  the  fracture  condition  on  the  hoop  stress  was  introduced, 
all  particles  have  a  negative  velocity,  i.e.,  toward  the 
center,  and  the  particles  at  the  7.5,  10,  and  12.5.  cm  stations 


29 


ARRIVAL  TIMES  AND  PULSE  WIDTHS  AS  FUNCTIONS  OF  RADIUS 


U 

w 

V) 

5. 


w 

w 


•J 


30 
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are  nearer  the  center  at  75  Msec  than  they  were  initially. 

In  many  of  the  runs  very  little  motion  has  occurred  at  the 
25  and  30  cm  stations  at  75  Msec.  Although  the  material  has 
attained  a  velocity  of  the  order  of  IQ3  cm/sec,  not  enough 
time  has  elapsed  for  any  significant  motion  to  occur. 
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4 .  SUMMARY 


Spherical  one-dimensional  computer  code  calculations 
have  been  performed  to  evaluate  the  importance  of  certain 
mechanical  material  property  parameters  on  stress  wave  propa¬ 
gation  in  dry  porous  rocks.  The  material  property  parameters 
were  chosen  in  a  range  that  is  representative  of  Nevada  Test 
Site  tuff.  It  was  found  that  for  a  given  boundary  condition 
the  attenuation  of  peak  stress  was  most  sensitive  to  the  ini¬ 
tial  porosity  of  the  material  with  the  attenuation  rate  in¬ 
creasing  with  increasing  porosity.  In  solid  material,  i.e., 
material  having  no  initial  porosity,  the  rate  of  peak  stress 
attenuation  increased  with  increasing  shear  modulus  and  with 
increasing  yield  strength  being  most  sensitive  to  the  shear 
modulus.  Shock  arrival  times  at  various  distances  from  the 
source  generally  increased  with  increasing  porosity,  thereby 
implying  that  average  shock  velocities  decreased  with  increasing 
porosity.  Particle  displacements  at  a  given  time  and  station  in¬ 
creased  with  increasing  porosity.  In  the  initially  solid  material, 
tensile  stresses  that  occurred  increased  with  increasing  yield 
strength  and  the  material  exhibited  in  some  cases  negative  dis¬ 
placements,  i.e.,  the  material  initially  moved  radially  outward 
from  the  source  and  then  pulled  back  inward  toward  the  source. 

This  pullrback  was  increased  by  subjecting  the  material  to  a 
fracture  condition  under  hoop  tensile  stress.  A  slight  decrease 
in  peak  stress  attenuation  in  porous  material  was  obtained  by 
increasing  the  value  of  the  hydrostatic  pressure  at  which  all 
voids  are  collapsed. 

The  work  reported  here  studied  the  effects  of  varying 
mechanical  material  property  parameters  in  a  dry  porous  rock. 
Further  parameter  studies  are  needed  using  different  equation- 
of-state  models  to  evaluate  the  effects  of  water  in  the  pores 
and  to  evaluate  thermodynamic  effects.  Parameter  studies 
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are  also  needed  on  hard  rocks  such  as  granit 
deformation  plasticity  theory  instead  of  the 
plasticity  based  on  the  infinitesimal  strain 


using  finite 
usual  theory  of 
approximation . 
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